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A large-scale cosmic magnetic field affects not only the growth of density perturbations, but 
also rotational instabilities and anisotropic deformation in the density distribution. We give a fully 
relativistic treatment of all these effects, incorporating the magneto-curvature coupling that arises 
in a relativistic approach. We show that this coupling produces a small enhancement of the growing 
mode on superhorizon scales. The magnetic field generates new nonadiabatic constant and decaying 
modes, as well as nonadiabatic corrections to the standard growing and decaying modes. Magnetized 
isocurvature perturbations are purely decaying on superhorizon scales. On subhorizon scales before 
recombination, magnetized density perturbations propagate as magneto-sonic waves, leading to a 
. small decrease in the spacing of acoustic peaks. Fluctuations in the field direction induce scale- 

' dependent vorticity, and generate precession in the rotational vector. On small scales, magnetized 

density vortices propagate as Alfven waves during the radiation era. After recombination, they decay 
■ slower than non-magnetized vortices. Magnetic fiuctuations are also an active source of anisotropic 

distortion in the density distribution. We derive the evolution equations for this distortion, and find 
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a particular growing solution. 
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O ' I. INTRODUCTION 

en 

. Recent observations reveal the widespread existence of magnetic fields in the universe and are producing much 
■ firmer estimates of their strengths in interstellar and intergalactic space. They also appear to be a common property 
^\ , of the intracluster medium of galaxy clusters, extending well beyond the core regions (see and references therein). 
0^ Strengths of ordered magnetic fields in the intracluster medium of cooling flow clusters exceed those typically associated 
with the interstellar medium of the Milky Way, suggesting that galaxy formation and even cluster dynamics are, at 
least in some cases, influenced by magnetic forces. Furthermore, reports of Faraday rotation associated with high 
redshift Lyman-a absorption systems seem to imply that dynamically significant magnetic fields may be present in 
condensations at high redshift [Q. The more we look for extragalactic magnetic fields, the more ubiquitous we find 
^ them to be. 

. Large-scale magnetic fields introduce new ingredients into the standard, but nevertheless uncertain, picture of the 
early universe. They seem unlikely to survive an epoch of inflation, but it is conceivable that large-scale fields and 
magnetic inhomogeneities could be generated at the end of that era or in subsequent phase transitions (see, e.g., [^). 
' Studies of magnetogenesis are partly motivated by the need to explain the origin of large-scale galactic fields. Typical 
spiral galaxies have magnetic fields of the order of a few coherent over the plane of their disc. Such fields could 
. . , arise from a relatively large primordial seed field, amplified by the collapse of the protogalaxy, or by a much weaker 
one that has been strengthened by the galactic dynamo. Provided that this mechanism is efficient, the seed can be 
as low as ^ lO^'^^G at present Q. However, in the absence of nonlinear amplification, seeds of the order of lO^^^G 
or even lO^'^G are required 

Determining whether the origin of galactic and cluster magnetic fields is primordial or post-recombination is a 
difficult task, since strong amplification in these virialized systems overwhelms any traces of their earlier history, fn 
contrast, magnetic effects on the cosmic microwave background (CMB) anisotropics, or any magnetic presence away 
from clusters and galaxies, can provide better insight into these early phases. If large-scale magnetic fields are present 
throughout the universe today, their structure and spectrum should bear clearer signatures of their past. Thus, 
improved direct observations, such as high resolution Faraday rotation maps and the study of extragalactic cosmic 
rays, may help in this respect For example, we would like to know whether or not the intergalactic voids are 
permeated by a widespread magnetic field, and whether there is magnetic field evolution in galaxies. If large-scale 
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magnetic fields were present in the early universe, were they dynamically significant, and if so, how have they affected 
the formation and evolution of the observed structure? It is known that element abundances constrain the strength 
of a primordial field at the nucleosynthesis epoch . Stronger limits on a primordial magnetic field are imposed via 
CMB anisotropics, since the field distorts the acoustic peaks and induces Faraday rotation in the polarization 

In this article we assume the existence of a large-scale ordered magnetic field of primordial origin a priori, and we 
investigate the magnetic effects on density inhomogeneities. Specifically, we analyze magnetized density perturba- 
tions, magnetized cosmic vortices (i.e., rotational instabilities), and magnetized shape distortion. Magnetized density 
perturbations were studied by Ruzmaikina and Ruzmaiki n [p| in Newtonian theory, while Wasserman fl^ looked at 
the rotational behavior of a magnetized fiuid. Kim et al. [|13[ derive a magnetized Jeans length, assuming that there 
are no density perturbations in the absence of the field. In a relativistic treatment, Battaner et al. investigated 
magnetized structure formation in the radiation era. Jedamzik et al. and Subramanian & Barrow have considered 
magnetic dissipative effects at recombination. 

We generalize aspects of these previous treatments by giving a fully relativistic analysis of the scalar and vector 
contributions of the magnetic field to the evolution of density inhomogeneity. We consider not only density pertur- 
bations and rotational instabilities, but also the shape-distortion effects of the field. Density perturbations are found 
explicitly in the radiation and dust eras, including a new solution that shows how the relativistic magneto-curvature 
coupling acts to enhance growth on superhorizon scales. The existence of the magneto-curvature coupling was first 
identified by Tsagas and Barrow | [T^ ]. The nonadiabatic magnetic effect on the modes is clearly identified, including 
the magnetized isocurvature modes. New solutions are also found for rotational instabilities, which are significantly 
affected by the field, and we show that magnetic effects actively generate shape distortion in the density distribution. 
We follow the relativistic analysis of cosmic electromagnetic fields given by Ellis [0, and we use the covariant and 
gauge-invariant approach to perturbations p"^-^2C|]. A covariant and gauge-invariant analysis of magnetized density 
perturbations was first developed by Tsagas and Barrow ^6 

We adopt the usual approximation that in the background, which is a spatially fiat Friedmann-Robertson- Walker 
(FRW) model,]] the ordered large-scale magnetic field is too weak to destroy spatial isotropy. The weak field approx- 
imation is an acceptable physical approximation when the field energy density is a small fraction of the isotropically 
distributed dominant energy density (see for further discussion). If one demands strict mathematical homogeneity 
in the background, i.e., if one refuses to accept a coherent test field in the background, then one must adopt a Bianchi 
I model for the background. However, this is a highly complicated approach, which in the end will give results that 
are practically indistinguishable from those with an FRW background. (We are currently completing calculations that 
confirm this statement (2^].) The standard assumption of very high conductivity is also made, so that we can ignore 
large-scale electric fields, while maintaining the desired coupling between the fiuid and the magnetic field. We use a 
single perfect fluid model, which is reasonable in the radiation era, but does not apply during recombination, while 
after last scattering, it means that our solutions only apply to a baryon-dominated universe. (See [2^ for a discussion 
of the effects of cold dark matter (CDM) potential wells on the field.) 

The assumption of a weak background magnetic field ensures that the field terms in our linearized equations are 
first-order. Current limits on the strength of a primordial magnetic field show that its influence is secondary relative 
to that of the dominant matter component. In practice, this allows for the possibility that second-order fluid terms 
can have a strength comparable to that of the linear magnetic terms. Nevertheless, such second-order terms do not 
contain any further information regarding the lowest order influence of the field on gravitational instability. We ignore 
the second-order fluid terms, even though they may be of comparable magnitude to the first-order magnetic effects. 
This approach is consistent at linear order, and allows us to isolate the lowest order magnetic effects on gravitational 
instability. Our aim is to identify the sources of the magnetic effects, calculate their impact to lowest order, and 
discuss their implications for the evolution of density inhomogeneities. 

In Sec. II, we outline the formalism and the main equations that govern the coupled evolution of density inhomo- 
geneity, the magnetic field and the curvature. Sec. Ill considers magnetized density perturbations, identifying the 
nonadiabatic effects of the field. We find a new solution on superhorizon scales in the radiation era, showing how the 
magneto-curvature coupling slightly enhances growth. In the radiation era, a small damping effect is wrongly pre- 
dicted when the magneto-curvature coupling is ignored. On subhorizon scales, magneto-sonic waves in the radiation 
era have a slightly increased frequency, leading to a decrease in the spacing of CMB acoustic peaks. In the dust era, 
the growing mode on small scales is slightly damped by magnetic effects. We also find the pure-magnetic density 
perturbations, i.e., the fluctuations created in a smooth plasma at magnetogenesis. These include growing modes. 
Magnetized isocurvature modes are characterized, and found explicitly on superhorizon scales. These modes all de- 
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cay, in the radiation and dust eras. Magnetized cosmic vortices are considered in Sec. IV. We show that magnetized 
vorticity is scale-dependent , and that the field generates precession in the rotational vector. We solve exactly for the 
rotational instabilities, showing that they propagate as Alfven (vector) waves on small scales during the radiation era. 
After recombination, such vortices persist for longer than non-magnetized vortices. Sec. V investigates magnetized 
shape distortion, showing that the field is as an active source of distortion. Purely magnetic distortion on superhorizon 
scales in the dust era is shown to have a growing mode. Conclusions are given in Sec. VI. 

We use units with c = 1 — SttG and our signature is ( — h H — h); spacetime indices are a,b, ■ ■ ■ and (square) round 
brackets enclosing indices denote (anti-)symmetrization. 



II. COSMIC MAGNETOHYDRODYNAMICS 



As noted above, the cosmic magnetic field B'^ must have weak energy density Pmag = ^BaB'^ to be consistent with 
observational limits, so that Ca ^ 1, where Ca is the Alfven speed. The Alfven speed, which effectively leads to a 
nonadiabatic increase in the sound speed, is given by 

4^^, (1) 

where p is the energy density of the cosmic fluid. In the limit of vanishing density inhomogeneity, i.e., in the 
background, the field is uniform, but its weak magnitude means that it does not disturb the background isotropy, so 
that the magnetic anisotropic stress is negligible in the background. In the actual inhomogeneous universe the field's 
influence propagates via: 

the background energy density and pressure (pmag = |/Omag) occur in terms of the form c^V, where 7-" is a 
perturbed quantity, and despite their weakness, they can have observable consequences (e.g., a change in the 
spacing of CMB acoustic peaks in the radiation era); 

spatial gradients of Pmag couple with gradients of p and thus alter the fluctuations of p (in particular, introducing 
nonadiabatic modes); 

the background direction of the field introduces anisotropy by picking out preferred directions in perturbed 
vector and tensor fields, and preferred directional derivatives of perturbed scalar/ vector/ tensor fields, leading 
to effects such as Faraday rotation; 

the background direction of the field is also the source of the magneto-curvature coupling, via terms of the form 
Kabcd.B'^, where Kabcd is the part of the curvature tensor which vanishes in the background; 

fluctuations in the direction of the field generate new anisotropics that can source magnetized vortices (leading 
in particular to Alfven waves) and shape distortion. 

We include all of these aspects in our analysis, so that we incorporate the full range of scalar (magnetic energy 
density and isotropic pressure) and vector (anisotropic stress) effects of the field, allowing for fluctuations in both the 
magnitude and direction. In order to provide a transparent relativistic generalization of Newtonian analysis, and to 
use variables that as far as possible have a direct physical interpretation, we adopt a covariant Lagrangian approach 
|p7|- p0| , p5 26 1 . This continues and develops the work of In particular, we discuss in detail the physical meaning 



and implications of the density perturbation solutions, and we extend the investigation to cover magnetized cosmic 
vortices and shape distortion. 

The cosmic perfect fluid defines a unique four-velocity Ua (with Uau"" — —1), and then hab — gab + UaUb, where 
gab is the spacetime metric, projects into the local rest spaces of comoving observers. The projection of a vector is 
^(a> = ha^Vb, and a projected second rank tensor Sab splits irreducibly as 

Sab = '^Shab + SabcS'^ + >5'(ah) j 

where S = habS""^ is the spatial trace, Sa = ■^SabcS'"^ is the spatial vector dual to the skew part of Sab, and 
S^ab) = [h{a'^hb)'^ — \h'^'^hab\Scd IS the projected symmetric tracefree (PSTF) part. Here Sabc = Vabcdu'^ is the 
projection of rjabcd, the spacetime alternating tensor. 

The covariant derivative splits into a comoving time derivative Ja - b = u'^^ cJa - bi and a covariant spatial derivative 
DcJa...& = h/ha'' ■ ■ ■ /ih^VdJe... ^ Then we define a covariant spatial divergence and curl that generalize the Newtonian 
operators to curved spacetime p5|: 
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CmWa = SabcD^'V" , CUrl5ah = ScdiaD'^Sb)''' ■ 

The fluid kinematics are described by the expansion Q = divw, four-acceleration Aa = Ua, vorticity uja — — |curlwa 
and shear aab = D(aUfc). Local curvature is described by the Ricci tensor Rab, while nonlocal tidal forces and 
gravitational radiation are described by the electric and the magnetic parts of the Weyl tensor, Eab — Ef^^ib) — 
Cacbdu'^u''' and Hab = H^ab) = ^£acdC'^'^beU'^ ■ The magnetized perfectly conducting fluid has energy density p and 
isotropic pressure p. The magnetic fleld is Ba ~ B(^a) i with energy density, isotropic pressure and anisotropic stress 
given respectively by 

Pmag — ^BaB'^ , Pmag — ^Pmag , T^ab — ~B^aBb) ■ (2) 

Then the total energy-momentum tensor is 

Tab = {P + Pmag) UaUb + {p + Pmag) hab + T^ab ■ (3) 

Notice that the absence of an electric field means that there is no energy flux (Poynting vector). The magnetic field 
appears from Eq. to behave like a radiation fiuid with anisotropic stress. However, this fluid picture does not 
fully encompass the vector properties of the fleld, and in particular, its coupling to the curvature. 

In the background, Ba is weak enough not to affect the isotropy, i.e. the anisotropic stress is negligible in the 
background, and Pmag ^ P- The background expansion is = 3iJ, where H = a /a is the Hubble rate. The 
background is covariantly characterized by 

D„e = r>aP = DaP = DfaSa = , (4) 
Aa=UJa=Q, (5) 
CTab = Eab = Hab = T^ab — T^ab — , (6) 

where 

Tlab = ha'^hb'^Rcd + Racbdu'^u'^ + DcUaDftU'' - QT)bUa , 

with Rab the Ricci tensor and Rabcd the Riemann tensor. Note that TZab is the intrinsic 3-Ricci tensor of spatial 
hypersurfaces only if tOa = 0; otherwise there are no such hypersurfaccs orthogonal to |l7|] . 

Quantities that vanish in the background are gauge- invariant, and they covariantly describe linear deviations from 
homogeneity and anisotropy. We collect below the linearized evolution and constraint equations given in ||l6|| , rewritten 
in the streamlined formalism of | p5[ |, which considerably simplifles the equations and facilitates analysis of their 
properties. The following covariant identities |Q are used in deriving the equations (assuming a flat background with 
vanishing cosmological constant): 

(7) 
(8) 
(9) 

(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

where the vectors and tensors vanish in the background, Sab — Si^ab) ■ The magnetic field itself does not vanish in 
the background, so that its projected derivatives do not commute at linear order: the vector identity in Eq. ( pT| ) is 
replaced by 

D[aDfc]-Bc = ^TZdcbaB"^ — Eabd^'^'Bc , 

where TZabcd is formed from Rabcd and the kinematic quantities p6[ |. This non-commutativity is the root of the 
magneto-curvature coupling found in n& . 



curl Da/ 


= -2fUa , 


(aBafy 


^aDaf + afAa, 


(D,/) 


= Ba (dV) + 2/curlu;, , 


{aDaJb...y 


= aBaJb- , 


B[a'Db]Vc 


-0-D[,Db]5^^ 


div curl V 


= 0, 


(div curl 5*)^ 


= icurl (div S')a , 


curl curl Va 


= D„(div1/)-DVa, 


curl curl Sab 


= lB^aidiyS)b)~B^Sab 
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A. Maxwell's equations 

In the infinite conductivity limit, Maxwell's equations [p^ , ^ provide three constraints, 

divB = 0, (16) 
cm\Ba^eatcB''A^+ja, (17) 
= i<7, (18) 

where ja is the current and q the charge density generated by fluctuations, and one propagation equation 

- -feSa + aabB'' + EabcB'^U^ , (19) 

which is the covariant form of the induction equation. Note that Bi^a) = Ba — Af,B''ua, and B"-Aa = to first order 
only in the case of a pressure- free perfect fluid |p^ . 

Contracting Eq. (^9|) with S", and neglecting the second order term cTafcTr"'', we deduce the radiation-like evolution 
law of the magnetic energy density, 

(B2)- + |e(B2) =0. (20) 
We can also derive the evolution of the anisotropic stress from Eq. (|9|): 

T^ab = —iHTTab - ^C^pCTab , (21) 

B. Conservation laws 

Energy density conservation is expressed via the equation of continuity 

p + e(l + u;)p = 0, 

where w = p/ p. Notice the absence of magnetic terms in this equation, since field energy conservation holds separately 
as a consequence of Maxwell's equations, as shown in Eq. (EOf). The two energy conservation equations imply 



[ciy = {w-\)Qcl. 

On the other hand, the field does enter conservation of momentum density 

(1 + w)pAa + c^D.p -I- EabcB^cmX B^^O, (22) 

where Cg = p/p is the adiabatic sound-speed squared (with D^p = CgDap). This equation reflects the momentum 
density exchange between the fluid and the fleld. The magnetic field is a source of acceleration (provided that curl 
Ba is not parallel to Ba); it can destroy the geodesic motion of the matter even in the absence of pressure. 

C. Kinematic equations 

Evolution of the expansion is governed by the Raychaudhuri equation 

e + 18^ + 1(1 + s.)p - ^^7^ + ^^^-1^ {2clA + clB) - A = , (23) 

where A is the cosmological constant, TZ = h°'^TZab is the projected curvature scalar, and 

A = aD''A., A. = ^, S.=D.S^ 
p B^ 

describe perturbations in the fluid and fleld energy densities. Note that the overall magnetic effect includes a coupling 
to the projected curvature, clTZ. 

Magnetic influence on cosmic rotation is encoded in the vorticity propagation equation 



5 



ioa + 2fl"wa = - jcurl Aa , (24) 
which may be rewritten, after ehminating the acceleration term via Eq. (|23) , as 



Ua + {2- icl) HoJa = -—^—B^'BbCmlBa . (25) 
^ ' 2(1 + w)p 

Thus there is a magneticahy induced vorticity component parallel to curl _Ba- The effect disappears if the directional 
derivative i3^D;,curl vanishes, i.e., when curlBa does not change along the magnetic force lines. 
Kinematic anisotropics evolve via the shear propagation equation 

&ab + 2fl"(T<jfc = — ^ rD(oAb) - — — — D(„Kt) + \'^{ab) 

a[l + w) ^ ' 2[l + w)p ^ ' 6(1 + w) ^ ' 

+ h-b + / , g"D,D^,Bb) - Eab ■ (26) 
(1 + w)p ^ ' 

The direct magnetic effects propagate through the field's anisotropic stress (TTab), as well as via anisotropics in the 
distribution of magnetic energy density iVii^a^b)) ^^^d of the field vector itself (D^^iSf,^). The latter effect vanishes when 
Y)(aBb) is invariant along the magnetic force lines. Also, the coupling between the field and the projected curvature 
has led to an extra magneto-geometrical contribution, c^TZi^ab)- 
The kinematic quantities also obey constraint equations: 

(diva)„ = |D,e + curios, , (27) 
divw^O, (28) 

Hab = CUrlcTah + D^aWf,) . (29) 



D. Curvature 

The electric and magnetic Weyl tensors obey Maxwell-like equations ]l7|,p6| 



Eab + SHEab - cml Hab = + w)aab + ^HWab , (30) 

Hab + SHHab + curl Eab = ^curl -Kab , (31) 

(div E)^ - iD,p +lBa-^ (div tt)^ , (32) 

(divH)^ ={l + w)puja, (33) 

where we have used Eq. (|2l|). For a magnetized fiuid the projected curvature tensor TZab is not in general symmetric, 
but has the form TZab = Tl{ab) + \Tihab + SabcTi-'^, where Q 

T^iab) = D(aA) - ^ {a^<Jab) + TTab , (34) 

TZa - curlAa - ^ (a^OJa)' , (35) 

7^ 2 (p- + A) (36) 



Note that the background relation 3H^ = p + A ensures that TZ vanishes in the background, so that it is gauge- 
invariant. By Eq. (M), the vector part of TZab is simply TZa = —Huja, which vanishes when the vorticity vanishes. 
The dimensionless curvature perturbation K, — a?TZ has comoving gradient ||l6|] 

oDalC = Ipa^I^a + o?Ba - ma^V)a^ , (37) 

which evolves as ||l^ 

{aDalCy = (2pc2Aa + aBa) + H^c^ cml lu a . (38) 

p(l + w) 
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E. Evolution of inhomogeneities 



The key gauge-invariant quantities describing inhomogeneity are the comoving spatial gradients of the fluid density, 
the field density, the expansion and the field vector: 

A„ = ^^^, Ba = IiaB\ Qa^aBaQ, Bab^aBbBa. (39) 

p 

Their propagation equations are 

A, = 3wHAa + il + w)ea + ^EabcB'^cm-l B^ , (40) 

P 

2 



-'^aeahcB^CMTlB" 



1 + w 



ai/curlwa , (41) 



Bab = -2HBab + Y"^j7^ [SB^a^b) + B^a^b]] - Bl^a^b) - -^[aQfc] 



+ aB^BcfTab + asab B^Bi^^Ud) ~ aB[aC\xv\uJb\ - aSacdB^H'^b ■ (42) 

Note how the magnetic field couples to the magnetic Weyl curvature via the last term in Eq . (p2[ ). 
By eliminating the expansion gradients from the time derivative of Eq. (^0|) , we arrive at |16| ] 

Ka = -{2 + icl - 6w) Hka + \[{}- 6c,^ + 8w - 3w^) p-2 (Sc^ - bw) A] A^ 

+ clB^Aa + f D^S, + - [{cl -w)p+{l + cl) A] ea,,S''curli?^ 
Ip P 

+ i(l + w)aBa + (^^y-) £a6cS''curlS<^> + [6(1 + w)c2 + 4c2] aiJcurlt^a • (43) 

In the Newtonian limit, Eq. ( ^ ) recovers the results of |ll]. The relativistic correction terms are the last three terms 
on the right hand side, i.e., the terms with Ba^ £abcB^cm:\I3'^'^\ and curltJa. 



III. MAGNETIZED DENSITY PERTURBATIONS 

The equations (^o|)-(^2|) provide the basis for a complete description of coupled density-magnetic inhomogeneities. 
We begin by isolating the evolution equations for the density perturbation scalars A (a covariant alternative to the 
density contrast 5p/ p) and B (describing fluctuations in the magnetic energy density), and the curvature perturbation 
JC: 

The required evolution equations are : 

A = - (2 + Sc^ - 6w) HA + \ [(1 - Qcl 8w - 3^^) p-2 (3c^ - 5u;) A] A 
+ clB^A icl + 2w)p-{l + icl) A] clB + \cIB^B 

+ \[{2- icl +?,w)p-{l + icl) A] clK. , (44) 

• 4 ; 4(d - w)H , , , 

■B=T77^ 7A + -^ ^A, 45 

3(l + w) H-w ^ ' 

JC^^^A + ^^B. (46) 
\ + w 1 + w 

This system of equations governs the coupling between density fluctuations A, magnetic fluctuations B and curvature 
fluctuations JC. Eq. (^6|) shows that K. grows if A and B are growing, while Eq. ( psf ) shows that if > w (which 
holds in the radiation and dust eras), then B grows in concert with growing A. 
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The magnetic field introduces a direct effect, via the term c^JC in Eq. (44), of the curvature on the density 
perturbations. In the non-magnetized case, there are two modes of A, which is governed by the single second-order 
equation (44) with Ca = 0. In this case, the evolution of A is independent of /C, and /C is determined once A is 
known, via Eq. (|4^). Magnetism introduces two additional modes, since the system has four degrees of freedom. 
These modes are nonadiabatic, and can source density perturbations, i.e., even when A(to) = = A(to), magnetic 
effects will lead to A 7^ for <> to- If one omits the magneto-curvature effect, then the evolution equation for /C, 
Eq. (^6|), is uncoupled from the system, which can then be decoupled via a third-order equation in A. Neglecting the 
magneto-curvature effect thus removes one of the additional nonadiabatic modes. 

For zero cosmological constant, we can solve the system analytically in the radiation and dust eras, treating super- 
and sub-horizon scales separately. Some solutions were given in There, however, magneto-curvature effects were 
neglected in three out of the four cases. Here, we generalize some of the solutions to incorporate the magneto-curvature 
coupling, and we show that the magneto-curvature coupling cannot in general be neglected, since it leads to important 
qualitative differences in the behavior of A. 



A. Radiation era 

During the radiation era, w = = ^, p = po{ao/a)'^, and the Alfven speed does not change along the fluid flow, 
i.e., Ca = 0, reflecting the radiation-like evolution of the magnetic energy density, as given by Eq. (|20|). For the Fourier 
modes with comoving wave- number k, we get 



A" = 



6' = A' , 



ao 



k 

kho 



-]lC' = A + ^ciB, 



A- 



k 

kho 



2 / \ 2 
a 

ao 



clB + 2clJC, 



(47) 

(48) 
(49) 



where a prime denotes d/ d{a/ao) and fcho = a^Ho is the comoving wavenumber of the horizon at oq. 



1. Superhorizon scales and the curvature coupling 



In the long wavelength limit k <C fcho, the system has the power-law solution 

a 



A = C, 



(0) 



JC 



E 



C(„) 



(50) 
(51) 
(52) 



where C^q) and C(a) are constants, and the parameter a satisfies the cubic equation 



a 



(53) 



The cubic has one positive and two negative roots. One of the negative roots corresponds to a decaying nonadiabatic 
mode. The other nonadiabatic mode is the C(o)-mode, which is constant. The remaining cubic roots correspond to 
the magnetized versions of the standard adiabatic modes, one growing and one decaying. Since is small, we can 
find the roots perturbatively. The zero-order roots are 0, —1, 2 (the a ~ solution is spurious in the non- magnetized 
case). To lowest order, we find that: 



o-cl + o{4) 

a={ -1 + 4 + 0(4) 



(54) 



•2 + 14 + (4) 
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Thus the adiabatic growing mode of the non-magnetized case is shghtly enhanced by magnetic effects (the enhancement 
is not felt to lowest order in c^); the adiabatic decaying mode decays less rapidly by virtue of magnetic effects; the 
decaying nonadiabatic mode decays very slowly; and the final, nonadiabatic, mode is constant. To lowest order 



C,i-) - 



-i+c: 



(0) 



C'(2-) 



(55) 



The magnetic and curvature fluctuations are given by equations (|51|) and (p2D, with a given by Eq. (|54|). 

This new solution in Eq. ( p5|) can be compared with the solution that arises when the magneto-curvature coupling 

\) falls away, leading to the quadratic — a — 



term c^/C is ignored in Eq. (|50|) [Iq ]. Then the last term in Eq. 
(2 - cl) = 0. To lowest order 



so that the density perturbation is given by 



(+) 



and the magnetic fluctuations are 



ao 



B 



-1 
2- 

+ Cr 



3^a ; 



3^a ) 



(-) 



ao 



3 ' 



(0) 



(0) • 



Clearly, omitting the magneto-curvature coupling has a significant qualitative impact. Not only is one of the 
nonadiabatic modes (C(2-)) removed, as expected, but we also find that the growing mode is slightly damped, at 
odds with the correct solution in Eq. (|55|). Thus the magneto-curvature coupling, which was identified in general in 
[ p^ , turns out to have a crucial role in increasing (even though it is only by a small amount) the standard adiabatic 
modes of density perturbations on large scales in the radiation era. ft is not reasonable to omit the magneto-curvature 
coupling in this case. 



2. Subhorizon scales and magneto-sonic waves 

At the opposite end of the wavelength spectrum, when k 3> fchO: we differentiate Eq. ( ^ ) and use Eq. (|4^) to 
decouple the system. Integrating once we get 



^',ao) ^ ^^(fcho) (ao 



; \ 2 / \ 2 

A; \ / a 



where Ck: is an additional constant associated with curvature effects. (We have ignored higher order terms in c^, 
given the weakness of the magnetic field.) This has solution (to lowest order in c^) 



A = 



where C(j) are constants. Si and Ci are the sine and cosine integral functions,^ and 



k a 
kho ao 
k a 
kho ao 



sin /3 



k a 
kho ao 



I k a \ 2 

cos 13- ) - Cbc^ 

kho ao , 



/3 = c, (1 + |c2) , 



(56) 



(57) 



Si(a;) = Jq t ^ sintdt and Ci(a;) — j + \nx + t ^(cost — l)dt, where 7 — 0.578 ■ ■ • is Euler's constant 
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where Cg = 1/ -v/S is the adiabatic sound speed. Thus (3 is the magnetized (nonadiabatic) sound speed of magneto-sonic 
waves. These waves differ shghtly in amphtude and frequency from the adiabatic acoustic waves. 

The magneto-curvature couphng, reflected in the nonadiabatic Cic mode, has the effect of shghtly modulating the 
amplitude of acoustic oscillations, with the effect decreasing as a/ao increases. The main magnetic effect is on the 



frequency. Comparing our result in Eq. (57) to the standard solutions of magnetic- free models (see, e.g., [UTi), we see 



that the field has increased the frequency of acoustic oscillations. Since a cx \/t, the magnetized acoustic frequency is 



Ho f ^ 



[l + ^cl) where = ^ I — ] . (58) 



This magnetic correction results from the the "tensioning" effect of magnetic force lines in the plasma, which produces 
a nonadiabatic increase of the sound speed via a contribution from the Alfven speed. As a result, the magnetic influence 
brings the acoustic peaks of short-wavelength radiation density oscillations closer, producing in principle an observable 
signature on CMB anisotropics An additional effect comes from the nonadiabatic constant mode in Eq. (^6|). Its 
presence suggests that the average value of the density contrast is nonzero, unlike the magnetic-free case. 



B. Dust era 



After recombination, in a baryon-dominated cold matter background, w = = c^, a = oq (^/^o)^/^ H = 2/3t and 
p — 4/3t^. The Alfven speed is no longer constant, but by Eq. (pG) varies as 



to 



2/3 



reflecting the fact that the magnetic energy density drops faster than that of nonrelativistic matter. Thus magnetic 
effects grow weaker as the expansion of the universe proceeds beyond recombination. 
The equations for the Fourier modes become 



A" 

B' 
IC' 



^) A' 



t 

i-Do 



^1 A 



2 l„2 
3 



8/3 



1 + i 



k 

kho 



3^ ■ 



^ (c^ 
3 V - 



5/3 



B, 



where a prime denotes d/d{t/to). Thus 



to 



2/3' 



B 



B^HA + Cb) 

^'-¥(^a)o(y)'''[A + CB] 



9 V a. 



8/3 



/c. 



(59) 
(60) 
(61) 

(62) 
(63) 



where = 0. We can now decouple the system. 



1. Superhorizon scales 

For long wavelength fluctuations, we get 

9 (^^-^ A'" + 36 (J-^ ' A" + 14 (^-^^ A' - 4A = , (64) 

to lowest order in c^. Note that curvature effects are quadratic in and do not contribute at this level. In fact 
) guarantee that, to lowest order in c^, curvature has no effect on magnetised disturbances in the 
dust distribution. We can solve Eq. (|6j), which is of Euler-type: 



equations (p9|), (61 
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2/3 ^ t \ ~'^ / 



A = C„(^-j +C,_,(^-J +C,„^^-j . (65) 

Thus, the field has simply added the nonadiabatic decaying mode C(2-) to the evolution of superhorizon density per- 
turbations, while the non-magnetized adiabatic modes are unchanged. The growth of large-scale matter aggregations 
proceeds virtually unaffected by the presence of the field or by curvature complexities. 

Magnetic effects on superhorizon scales in the dust era do not change the adiabatic growing mode to lowest order 
in c^. The adiabatic decaying mode is also unchanged (unlike the radiation case). However, a new nonadiabatic 
decaying mode arises, which decays less rapidly than the adiabatic mode. 



2. Subhorizon scales 



On subhorizon scales 



A' 



1 - ^ icl^ 



9 V'-a/o I I, 

\ "-hO 

where again we have ignored terms of higher order in (c^)o. The solution is 



A = C, 



(+) 



C(-) 



C(H-) 



to 



-2/3 



4(c2)oF-9fc2, 



with 



a± 



-l±5Wl-f|(c2^ 



a/0 



(66) 



(67) 



The magnetic infiuence is expressed in two ways: additional decaying (C(e_)) and constant (Cg) nonadiabatic modes; 
and modification of the non- magnetized adiabatic modes (C(±)). The net effect is to inhibit the growth of matter 
aggregations, as noted also in the Newtonian case Note that the magnetic effects, direct or indirect, become 

less important after matter-radiation equality, due to the decrease of the Alfven speed. The damping of the growing 
mode is greater on smaller scales. Indeed there is a minimum scale, below which the solution in Eq. (p^ ) oscillates, 
since the magnetic pressure balances gravitational infall. This magnetic Jeans scale follows from Eq. (|67[): 



Am J (to) = |7r\/6Aa(<o) , 



(68) 



where 



Cat — gCaAh 



is the Alfven horizon, with Ah = H^^ the Hubble scale. In fact, given the weakness of the magnetic field, it is likely 
that kinetic pressure cannot be ignored near the magnetic Jeans scale. In this case, a more sophisticated analysis is 
necessary, to incorporate nonrelativistic pressure effects in baryonic matter. On scales well above the Alfven horizon 
(equivalently, magnetic Jeans scale) but well within the Hubble horizon, i.e., for 



fcho < fc < /cao where fca = -z—k\, 

2Ca 

we find that the magnetized corrections a± of the adiabatic exponents are 
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The way in which magnetic effects act to increase the adiabatic Jeans length may be quahtatively understood as 
follows. Consider a tube of magnetic force-lines with instantaneous cross-sectional area SS. In a perfectly conducting 
medium the field remains frozen into the fluid, i.e., the magnetic force-lines always connect the same particles [ pT| . 
More precisely, the induction equation (|l^) shows that a^B°- is a connecting vector. Thus the volume of the tube is 
given by 

6V = d£6S cx a^BSS. 
However, we also have that in general, 6V oc a^. It follows that 

{BSSy =0. 

The conservation law in Eq. ( pO[ ) shows that B (x a^^; thus 

6S (X , 

so that the cross section of the flux tube increases as the expansion redshifts the energy density of the field. Thus the 
field acts against gravitational infall. 



C. Pure- magnetic and magnetized isocurvature perturbations 

We have seen that the magnetic field introduces nonadiabatic modes in the density perturbations. This means that 
the field itself can generate fluctuations in the density, even when there are no primordial density fluctuations. Thus, 
if 

A{to) = = A{to) , (69) 

where Iq is the epoch of magnetogenesis in the early radiation era, then nonzero A will arise purely from the magnetic 
field; in the absence of magnetogenesis, Eq. ( |69|) would imply A — for t > to- These nonadiabatic pure magnetic 
density perturbations can be found explicitly from the solutions given above. On superhorizon scales (assuming that 
the field is created on these scales at to), the general solution in Eq. (p5[ ) implies with the initial conditions in Eq. 
(|69|) that the pure-magnetic nonadiabatic mode is (to lowest order in c^) 



-3+c„ 



+ C(o)+C(2_)(^^j . (70) 

The pure-magnetic density perturbations have a dominant growing mode of the same strength as in the non-magnetized 
adiabatic case. The decaying modes are in fact the isocurvature part of the pure-magnetic density perturbations, as 
we now show. 

Equation (|6^) is often taken to characterize isocurvature perturbations, but it does so only in specific cases [ pO[ . 
For magnetized perturbations, this is not the isocurvature condition. Isocurvature density perturbations are those 
for which the curvature perturbation of the initial hypersurface orthogonal to the fluid flow is spatially constant, i.e., 
{aDal^){to) — 0. (There is also the implicit condition that uja — 0, which is necessary for the existence of the spatial 
hypersurface.) Taking the the comoving divergence of Eqs. (|3^) and (pO|), we find the condition for magnetized 
isocurvature perturbations: 

A+^l-w)HA^l{l-w)clHB-clHJC at t = to . (71) 

In the non-magnetized case Ca — 0, it is clear that this condition is satisfied by Eq. (|69|), but when Ca > 0, then Eq. 
( |69| ) does not characterize isocurvature perturbations. 

As an example, consider the implication of the magnetized isocurvature condition in the dust era, on scales well 
above the Alfven horizon but well within the Hubble horizon, i.e., k^o <C fc <C fcao- Then Eqs. ( [7l| ) and (|6|) give, to 
lowest order, 



C(+) = f (ca)o [C(-) + C'e] - ^Cfc 
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On superhorizon scales, Eq. ( |7l| ) holds for all t by virtue of Eq. (^8|), which implies (aDa/C)' = 0. The magnetized 
isocurvature condition then selects a sub-class of the general superhorizon solutions found above. In the radiation 
era, we find (to lowest order in c^) 



C(2-) I — 



to 



Equation ([7 
dust era, 



(72) 



arises from the general superhorizon solution Eq. (|55D by eliminating the non-decaying modes. In the 



C(2-1 I — 



to 



(73) 



which is a special case of the general superhorizon dust solution Eq. (pq), once again without the growing mode. 

It is clear from Eqs. (|7^ ) and (|7^) that magnetized isocurvature perturbations on superhorizon scales are purely 
decaying. They are very different from the pure-magnetic solution Eq. (^o|) that is based on the initial conditions in 
Eq. (69). The latter has a constant and a growing mode. Magnetized nonadiabatic perturbations on superhorizon 
scales can contribute to the growing mode and generate a constant mode, whereas the magnetized isocurvature 
perturbations are purely decaying. 



IV. MAGNETIZED COSMIC VORTICES AND ALFVEN WAVES 

In the previous section, we generalized the results given in [Q, which itself provided a relativistic extension of 
previous work on magnetized density perturbations. A general inhomogeneous perturbation is characterized not only 
by its magnitude, i.e. the density perturbation A, but also by its rotation and deformation properties, as described 



in general terms in m9 20 |. Recently, these properties were investigated in CDM, and it was shown how the small 



stresses (isotropic and anisotropic) from residual velocity dispersion can have an important effect on rotation and 



deformation, even though the effect on density perturbations is effectively negligible 1 29 

The evolution equations for rotational and deformation variables in an imperfect fluid were derived in ]20( |. The 
evolution equations for inhomogeneities were coupled to causal transport equations for viscosity and heat conduction. 
By Eq. (^), a magnetized perfect fluid can be considered as an imperfect fluid with anisotropic stress, and the 
equations of may be specialized to this case. However, the system needs to be completed by evolution equations 
for the magnetic stress, which are determined by Maxwell's equations. Here we investigate the coupled equations 
governing rotational and deformational inhomogeneity in the fluid and magnetic fleld. 
The comoving gradient of the density inhomogeneity splits irreducibly as 

aBbAa = ^Ahab + SabcW + ^ab ■ 

The density perturbation is the comoving divergence, the rotational part is given by the comoving curl and the 
deformation part is the comoving PSTF derivative: 

A = aD'^Aa , Wa - -iacurl Aa , Cah = aD(^Ab) . 

The vector Wa governs rotational instabilities in the density distribution of the matter, and div = 0. On the other 
hand, £,ab determines the volume-true anisotropic distortion, with h°'^^ab ~ 0. Both quantities describe differential, 
i.e., infinitesimal, properties. Here we focus on rotation, and in the next section we look at anisotropic deformations. 

A fundamental property of rotational perturbations is that they are proportional to the vorticity vector. This arises 
from the identity Eq. (^ which ensures that the curl of any gradient field, such as A^, derives from the vorticity. It 
follows that 

Wa = -Sa^Hil + w)u;a ■ (74) 

Not only Wa, but also rotational perturbations in magnetic density and expansion inhomogeneities, are parallel to 
the vorticity: 

, 4 H 

cuilBa = curlD^S^ = — — -Wa , curlOa = curlaD^e = -— rrrWa ■ 

3(1 + w) (1 + w)H 
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Thus all rotational instability in the magnetized medium arises from vorticity. As we have seen from the vorticity 
propagation equation (|2^), magnetic inhomogeneities can source vorticity. Rewriting Eq. (^ ) using Eq. (Q), we 
have 



B''Dbcm-\Ba . 



(75) 



Thus the field is a source of vorticity provided that its curl varies along its force lines. Furthermore, the effect of 
the field is to induce precession of the rotational vector Wa- In the absence of the field, Wa remains parallel to Wa, 
so that the initial direction is preserved along the fluid flow. By contrast, in the magnetized case, Wa is no longer 
parallel to Wa, and the initial direction changes along the fluid flow. The rate of precession is 



■^prcc 



iB^DtcmlBal 

2{l+w)p\LJa\ 



(76) 



Equation ([75|) shows how the magnetic field can become a source of density vortices. However, the field effect upon 
pre-existing rotational perturbations is not clear yet. To quantify the magnetic influence on Wa we need to go one 
step further and obtain a decoupled equation for the evolution of Wa ■ We take the curl of Eq. ( p3| ) , using the above 
results and the identities in Eqs. (^, (10|) and (|lj), and we arrive at the required evolution equation (with A = 0): 



Wa + {4.- 3w) HWa + l[l~7w + 3c2(l + w)] pWa 
This is a wave equation for Wa, with signal speed Va given byQ 



3(1 + w) 



(77) 



3(1 + w) 



Propagating solutions of this equation are Alfvcn waves (compare |^), i.e., incompressible, vector waves, as opposed 
to the compressible, scalar magneto-sonic waves. In the non-magnetized case, the signal speed vanishes, and no wave 
solutions exist. We note also that the only scale-dependence arising in the wave equation is via the magnetic term. 
Thus magnetized vortices are scale-dependent, unlike the non-magnctizcd case. 
Decomposing the solenoidal vector Wa into Fourier modes W, Eq. ( |77| ) gives 



W + (4 - 3w) HW + U [1 - 7w + 3c2(l + w)]p + 



3(l + u;)a2 



W = 0. 



(78) 



Clearly, Wa can only grow if the term in square brackets becomes negative. If zi; = 0, as in the radiation and dust 
eras, then the quantity 1 — 7w -I- 3(1 -I- w)cl = {Sw — l)(w — 1) becomes negative when ^ < w < 1. Thus it is only 
when matter stiffer than radiation dominates the universe (and is coupled to the magnetic field), that vortices in the 
density distribution can grow (in the linear regime). Furthermore, the presence of the magnetic field ensures that 
such growth occurs only on scales larger than the critical "rotational Jeans" wavelength 



A 



rj 



27rCa 



3(3w;-l)(l-w2)p 



1/2 



which is small due to the weakness of the field. 
In the radiation era, Eq. (uSt) becomes 



W = 0, 



where a prime denotes d/d{a/ao), and fca 
the general solution 



2k\i/ca is the wavenumber of the Alfven horizon Aa = ^CaAh. This has 



^ A similar equation was derived in for a fluid with shear viscosity; in that case = ri/[rp{l + w)], where r] is the viscosity 
and T is the causal relaxation time. 
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which describes Alfven waves. The Alfven frequency 



fcaO 



(80) 



where Jz/ac — i^ac.mag — t'ac IS the excess magnetic acoustic frequency given in Eq. ([S^)- Thus local differential vortices 
in the density distribution are "flip-flopping" in concert with the acoustic oscillations in the density perturbations. 
Alfven waves are a purely magnetic effect, arising from the fluctuations in the magnetic fleld direction. These waves 
have decaying amplitude, in common with non-magnetized (and non-propagating) vector perturbations. 

On scales well beyond the Alfven horizon, k <C fcao, the oscillatory behavior is not felt, and Eq. ( |7^ ) gives, to lowest 
order, 



W = C(i) 



oo 
a 



k 



(2) 



k 
fcaO 



On superhorizon scales, the oscillations disappear: W — > C(i)(ao/a), regaining the standard non-magnetized result. 
Thus, before matter-radiation equality, and on scales much larger than the Alfven horizon, density vortices evolve 
unaffected by the presence of a cosmological magnetic field. 

After equality, in the matter-dominated dust era, Eq. (|7^) becomes 



3(^1 W" 



where a prime denotes d/d{t/to). The solution is 



k 

fcaO 



W = 0, 



(81) 



where 



a± 



-5± Wl - 12 



(82) 



(83) 



Therefore, any rotational instabilities present in the density distribution of the dust die away with time, as they do in 
non-magnetized cosmologies. The field effect on a given mode k is to reduce the depletion rate of Wa by an amount 
proportional to the initial Alfven speed squared, (0^)0- Thus, magnetized dust universes will contain more residual 
vortices than magnetic-free ones. However, the effect is confined within a narrow wavelength band beyond the Alfven 
horizon Aa. On much larger scales, the field influence becomes negligible, and Wa oc to/t as in non-magnetized models. 
On scales with k > fcao/'v/TS, i-e. within a few times the Alfven scale, Eq. (|83|) shows that the density vortices oscillate 
as Alfven waves. 



V. MAGNETIZED SHAPE-DISTORTION 



We monitor anisotropic deformation (shape distortion) in the density distribution of the medium through the PSTF 
tensor ^ab = aD(a^6)- This is associated with density variations that do not represent matter aggregations, since 
the associated divergence of Aq is zero, but rather describe changes in the local anisotropy pattern of the density 
gradients. 

Distortion in the density is coupled to distortion in the expansion and the magnetic energy density, defined via the 
PSTF tensors 

dab^a^^(a^^)Q, P^, = ^B^aT)h)B\ (84) 

which vanish in the background and are thus gauge-invariant. The propagation equations for ^at, ^^ab and Pab follow 
from the comoving PSTF derivatives of Eqs. (|4C|)-(^): 



15 



Cafc = SwH^ab - (1 + W)^ab + jciH(3ab " cfi/K^b - 3ff/^a6 , (85) 

1 + w 2(1 + w) 



4c2 

^clpKab - ^PPab ~ ( 6c^ + j ^^^^ab , (86) 

• 4 • 4(c2-u;)i7 

Pab = xt:— r?ab H — Ub ■ (87) 

3(1 + w) 1 + w 

The additional gauge-invariant PSTF tensors 

Kab = a^'Tl{ab) , '^ab = a^D(QCUrltJf,) , Pab = -B^DcBi^ab) 1 

respectively describe distortions caused by projected curvature, rotation and by anisotropics in the distribution of the 
magnetic field gradients. The first is due to the natural coupling of the field to the curvature and is given by 

Hab = (^TTafc — Hdab + Eab) , 



obtained from Eq. (|3J) by means of the shear propagation equation (^6|) . The second arises from the fluid flow, which 
generally is not hypersurface orthogonal. It has no impact on deformation if the rather special condition D^^iCurlw^) = 
holds. Finally, the effect of pab vanishes when any anisotropics present in the distribution of Bab = ciDtBa remain 
invariant along the magnetic force-lines, that is when B^T)cBi^ab) — 0. Note that both the scalar and vector aspects 
of the field contribute to shape distortion. 

Equations (|8^) and (^6|) combine to provide a second order differential equation, also obtained by taking the 
comoving PSTF derivative of Eq. (Esf). With A = 0, we have 



Lb - - (2 + icl - &w) Hiab + Hi - + 8w - 3w') pU + cl^'^U 

-i (1 - + 2w) clpPab + ^clD^f3ab + 5(2 - + 3u;) {clpKab + Sppab) 
+ [6(1 + wy, + 2cl] HWab- 

This is coupled to Eq. ( p7| ) for the growth of infinitesimal distortions in the magnetic energy density. The other 
source terms, namely Kab, jJ'ab and Wab, evolve according to 

1 • {c^ + 3w)H^ '^clH ^ AH ^ ^2 ,on^ 

l^ab = —— S.ab H — S,ab + — Pab — — /^ab " ^TUab — F) (Taf, , (89) 

I + W I + W 1+W 1 + W 



3(1 + w)^"'' 1 + w 

ctH ^ ctH 



flab = - [4 - 3(1 + W)] Hpab + " M b + \^ U 

Sil + w) 1 + w 



2(l + w)^-+3(rf^^- + ^^^«'°''^-' (90) 

Wab = -(2 - 3c^,)Hmab + — ^— DeD(,i?'^D|,| [Bb)^ - B',^] . (91) 

2(1 -I- w)p ^ I I L / /J 

We used the propagation equations (pw and (E^) for Eab and iTab, which imply 



Eab = -SHEab + ^H-Kab " ^(1 + w) pO ab ~ D^(T<jb + — {."^ ab + '2vJ ab) , 



on using the constraint equations (|^ and (p9|). Equation ( |90| ) requires Eqs. ( p^ ) and (|8^). Finally, to obtain the 
evolution formula of Wab we have successively taken the comoving curl and the comoving PSTF derivative of Eq. ([2^). 

The system of equations (^7[)-(|9l|) provides in principle a complete description of linear infinitesimal shape distortion 
generated by magnetic effects, provided we have a prescription for the Ti^Uab terms in Eqs. (|9|) and (|o|), and for the 
last term on the right of Eq. (^). Even without these terms, the system is too complicated to analyze in general. 
However, it is clear in general terms how magnetic effects will actively generate distortion. We can illustrate this by 
comparing with the non-magnetized case in a simple example. 

For simplicity, consider superhorizon scales in the dust era (neglecting vorticity). Suppose that at a given event 
(to, xo), we have no initial distortion or rate of distortion: 
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iU)o = = (Cafc)o , (92) 
In the non-magnetized case, the distortion system coUapses to the single equation 

tab — —'^H^ab + ^H'^^ab ■ 

It follows that along the fluid flow line through xo, no distortion is generated: 

S.ab{t,Xo) = 0. 

The evolution is purely passive^ or inertial, i.e., distortion can only develop if it is there a priori. In the magnetized 
case, by contrast, distortion is actively and nonadiabatically generated by magnetic effects. Equation (88) shows that 

Ho'^{iab)o = (Ca)o [-| (Kafc)o + 2 (Kah)o] + 6 {^lab)a , 

so that {^ab)a 7^ 0. Distortion is immediately generated along the flow line. In fact, the distortion has a growing 
mode, as we now show. 

The superhorizon scalar modes of the distortion tensors satisfy a system that follows from Eqs. (^ ) and (^8|)-(^0|): 



a' 



4 H 



3 V-^a^o 



5/3 



9 V-a 



8/3 



[3/3 - 4k] 



2/3 



to 



(93) 
(94) 

(95) 
(96) 



Equations (M) and (Pq) integrate to 



2/3 



(97) 



where = = F^. Then Eq. 



transforms Eq. ( |95|) into 



to 



2/3 



(98) 



According to Eq. (p7|), the effect of any anisotropics present in the distribution of Bab = aDaBt on decreases after 
matter-radiation equilibrium. Since /i is a key source of shape-distortion, we expect the evolution of ^ to approach 
that of A as the universe expands. Equation (£3) gives 



9 - C'" + 36 - e" + 14 - k' - 4^ = 



on using Eqs. (|97|) and (|98|). This has the same form as the corresponding density perturbation equation (|6J), and 
the solution is thus of the form Eq. (|65|). Imposing the initial conditions in Eq. (92), we find that 



to 



2/3 



to 



-2/3' 



(99) 



where F is a constant. Thus the shape distortion has a growing mode along the fluid flow line, due purely to magnetic 
effects; in the absence of the magnetic field, F = 0. (A similar situation arises in the simpler case of distortion 
generated by velocity dispersion [p9[.) 
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VI. CONCLUSION 



We have given a fully general relativistic treatment of the scalar and vector effects of a weak large-scale magnetic 
field on cosmological density inhomogeneity. This refines the results of on magnetized density perturbations, and 
extends that work to analyze magnetized vortices and shape distortion in the density distribution. Our covariant 
Lagrangian approach allows us to derive gauge-invariant evolution equations for all these aspects of density inhomo- 
geneity in the general linear case, i.e., incorporating all fluctuations of, and couplings between, the field, the fluid and 
the curvature. In summary, r nag netized density perturbations are governed by Eqs. ( p3)- ( |4(:), magnetized density 
vortices are governed by Eq. (JTfj), and magnetized shape distortion is governed by Eqs. (|87|)~(3l|). 

We give the solutions in closed form for magnetized density perturbations and vortices, in the radiation and dust 
eras. Some of the scalar solutions and all of the vector solutions are new. For magnetized shape distortion, we found 
a special solution with a growing mode. 

Given the overall weakness of the field, the magnetic effects described here are secondary relative to those of 
the matter. In some cases, second-order fluid effects may be comparable in strength to first-order magnetic effects. 
However, their presence does not affect the field impact, which remains the same. Thus we can neglect the second- 
order effects in a consistent linear analysis that probes the lowest order magnetic effect on density inhomogeneity. To 
lowest order, the magnetic influence on gravitational instability may be summarised as follows: 



1. The magneto-curvature coupling, which is a direct consequence of the fleld's vectorial nature, first identified in 
[ p^ , has an important influence. In particular, on superhorizon scales in the radiation era, this coupling slightly 
enhances the growing mode of density perturbations, as shown in Eqs. (|5^) and (^J): 



ao 



2 '-a 



C(i-) - 



+ C'(o) + G(2-) — 



When the coupling is neglected, the growing mode is incorrectly found to be damped relative to the non- 
magnetized case. 

2. Magneto-sonic waves in the radiation era are given in exact form in Eq. (|5^), 



A 



. ( k a 
sm /3- 

cos [(3^— ] - CbcI 
\ Kho ao , 



This shows the nonadiabatic modulation of the amplitude and increase in the frequency of acoustic oscillations. 
These effects, together with the nonzero average value implied by the Cb term, have potentially important 
implications for the CMB acoustic peaks, some of which have been investigated in 

3. In the dust era, subhorizon magnetized density perturbations are given exactly in Eq. (|6^), which leads to the 
magnetized Jeans scale in Eq. (|68|): 

Ainj(io) = ^7r\/6(ca)o Aho ■ 



On scales such that Amj ^ A ^ Ah, the density perturbations are 



where e — |(/c/A:ao)^- This shows the small damping effect on the adiabatic growing mode, as well as the 
new nonadiabatic modes. These results imply small modifications to structure formation in the linear regime. 
However, they are limited by the fact that we have neglected any non-baryonic matter or cosmological term. 

4. Pure-magnetic density fluctuations, which are induced in an initially smooth fluid by magnetogenesis, are given 
on superhorizon scales by Eq. (|70|). This solution would be important in any attempt to model large-scale 
structure formation as seeded by magnetogenesis. 
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5. Magnetized isocurvature perturbations are characterized by Eq. (|7 
purely decaying. 



On superhorizon scales, these modes are 



6. The field is a source of incompressible rotational instabilities, and the condition for this to happen is given via 
Eq. ([75]). Magnetized density vortices are shown be scale-dependent and to precess, at a rate given by Eq. j?^). 

The general propagation equation for these vortices (i.e., incorporating all relevant effects) is given by Eq. (|77 



Wa + {4.- 3w) HWa + i [1 - 7w + 3c2(l + w)] pWa = 



3(1 



In the radiation era, the Alfven wave solutions are given exactly in Eq. (^ 
wave-speed are 



The Alfven frequency and 



Ho 

TT 



These results generalize some of the theoretical results of Durrer et al. who then go further and apply the 
results to determine the effect of Alfven wave modes on CMB anisotropics. 



7. After recombination, magnetized density vortices are given exactly in Eq. (82). They decay like their adiabatic 
counterparts, but at a slower rate, so that rotational instability persists for longer in a magnetic universe. This 
will have a small effect on structure formation in the linear regime. 

8. Finally, we have investigated for the first time magnetic effects on infinitesimal shape distortion in the density 
distribution. The magnetic influence is manifold. Anisotropics in the field energy density, together with those 
in the distribution of the magnetic vector itself are direct sources of density deformation. The field's coupling 
to curvature and rotation also acts as an indirect source of magnetically induced shape distortions. Following 
the evolution of shape-distortion along the worldline of a fluid element, we showed that the field is an active 
source of distortion. On superhorizon scales, we showed via a special solution of the shape-distortion system 
that there is a growing mode of magnetized shape-distortion [see Eq. (p9|)]. 

Unlike the magnetic effects on density and rotational perturbations, which are small corrections of the non- 
magnetized results, magnetic effects on shape-distortion constitute a significant change from the non-magnetized 
(and passive) case. (A similar statement applies in the case of velocity dispersion effects in CDM ^^.) The 
results on magnetized shape-distortion have potentially important implications for (linear) structure formation. 
Not only is distortion actively generated once scales re-enter the Hubble horizon and begin to collapse, but it 
is also actively generated while the scales are beyond the horizon. Of course, the shape distortion in the linear 
regime will be overwhelmed by effects that arise during the nonlinear stages of collapse. 
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